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By growing monolayer FeSe on SrTiO 3 (001) surface, researchers obtain the highest superconducting 
transition-temperature for iron-based superconductor. Here, we study the antiferromagnetic (AFM) 
checkerboard monolayer FeSe adsorbed on SrTiO 3 (001) surface. We show that the system has a considerable 
charge transfer from SrTiO 3 (001) substrate to FeSe monolayer, and so has a self-constructed electric field. 
The FeSe monolayer band structure near the Brillouin zone (BZ) center is sensitive to charge doping, and the 
spin-resolved energy bands at BZ corner are distorted to be flattened by the perpendicular electric field. It is 
revealed that the striking disappearance of the Fermi surface around the BZ center can be well explained by 
the AFM checkerboard phase with charge doping and electric field effects. We propose a tight-binding 
model Hamiltonian to take these key factors into account. We also show that this composite structure is an 
ideal electron-hole bilayer system, with electrons and holes respectively formed in FeSe monolayer and Ti0 2 
surface layer. 

Due to its simple crystal structure, prominent antiferromagnetic magnetic ordering, and significant pressure 
effect on superconductivity, during the last few years FeSe binary superconductor has been drawing 
extensive attention from fundamental mechanism of Febased superconductivity to promising application 
in such areas as superconducting wires and thin films 1 . Recently, monolayer FeSe with atomic flatness was 
successfully grown on SrTiO 3 (001) surface 2 , and strikingly, the high-temperature superconductivity signature 
from the monolayer was revealed by high energy resolution scanning tunneling spectroscopy (STS). It is even 
hopeful that with further improving the sample quality, the transition temperature is probably as high as 77 K, the 
liquid nitrogen boiling point. The subsequent angle- resolved photoemission spectroscopy (ARPES) experiment 3 
shows that the Fermi surface of this FeSe monolayer only consists of electron pockets near the Brillouin zone (BZ) 
corner, with no hole pockets around the zone center observed. The T c = 65 ± 5 K is realized by optimizing the 
annealing process 4 , and spin density wave is also identified in FeSe/SrTi0 3 thin films 5 . Theoretical efforts are 
being paid to explain these intriguing STS and ARPES reports. The soft phonon of SrTiO 3 (001) was found to 
strengthen the Cooper pairing 6 . S 4 symmetry was also noticed 7 , wherein it was suggested that cuprates and iron- 
based superconductors share an identical high-T c superconducting mechanism. The density functional theory 
(DFT) total-energy calculations gave that collinear antiferromagnetic order (CAFM) FeSe monolayer on Ti0 2 
terminated SrTiO 3 (001) surface is the most stable structure, and there is neither hybridization nor charge transfer 
between FeSe monolayer and SrTiO 3 (001) surface 8 . The calculated Fermi surface for this CAFM structure exists 
around the M point in the folded BZ, and there is a Dirac-cone-like bands of FeSe near Fermi energy, which 
unfortunately are not observed by the ARPES measurement 3 that the Fermi surface disappears around the zone 
center. A more recent DFT study 9 on a single FeSe monolayer (without substrate) showed that an AFM checker- 
board phase (we will use AFM for simplicity in the following texts) instead of CAFM phase can have the Fermi 
surface resembling the ARPES data. Basically, however, the interaction between AFM-ordered FeSe monolayer 
and SrTiO 3 (001) surface remains largely unclear in literature, despite its urgent importance on understanding the 
key role played by the interface in controlling the Fermi surface topology and high-temperature superconduct- 
ivity of FeSe monolayer. 

In this work, we carefully investigate with DFT based ab initio calculations the electronic structure of AFM FeSe 
monolayer on Ti0 2 terminated SrTiO 3 (001) surface. We show that the calculated electronic structure of the 
system is close to the ARPES experimental results. The charge transfer and electric field are induced by the 
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Figure 1 | The side view (a) and top view (b) for the atomic structure of 
monolayer FeSe absorbed on Ti0 2 terminated SrTiO 3 (001) surface. The 

magnetic ordering is also represented in top view (b) as red arrows. 

interaction between SrTiO 3 (001) substrate and AFM FeSe mono- 
layer, and turn back out to play an important role in modulating 
the band structure and Fermi- surface topology of the FeSe mono- 
layer. One extraordinary consequence is that the spin-resolved elec- 
tron pockets surrounding the BZ corner are distorted and become 
flattened. We propose a tight-binding model Hamiltonian that can 
describe these key factors. 

Results 

We use 6-layer SrTiO 3 (001) slabs to mimic the surface of SrTi0 3 
substrate, as shown in Fig. 1(a). The SrTiO 3 (001) surface is either SrO 
or Ti0 2 terminated with nearly the same surface energy 10 " 12 . 
According to the experiment 2 , we use the Ti0 2 terminated surface 
in our atomic models. FeSe monolayer is adsorbed on the 1X1 two- 
dimensional unit cell of Ti0 2 -terminated SrTiO 3 (001) surface. To 
test the influence of SrTi0 3 antiferrodistortion on surface electronic 
structures, we also calculated a \fl x y/l supercell, in which the lower 
four layers of Ti0 2 and SrO are fixed to their bulk positions and the 
octahedral rotation angle is set to the experimental value 2.1° 13 . In 
our models, the vacuum space is larger than 10 A. The results are also 
checked for vacuum space larger than 20 A. The bulk calculation 
shows that the lattice parameter of SrTi0 3 is 3.901 A, which agrees 
well with experiments and previous theoretical results. The adsorp- 
tion structure relaxation is performed with fixed lattice parameter of 
3.901 A. 

The relaxed adsorption structure of FeSe monolayer on Ti0 2 ter- 
minated SrTiO 3 (001) surface is shown in Fig. 1. The lower layer Se 
atoms are on the top of Ti atoms, while the Fe atoms locate on the top 
of O atoms. The vertical distance between Ti and Se atoms is 3.13 A, 
which is slightly larger than that of the CAFM phase by 0.07 A. The 
vertical distance between Fe and O atoms is 4.43 A. The O atoms in 
the Ti0 2 top-layer arise slightly, they are 0.1 A higher than Ti atoms. 
The distance between top (bottom) Se atoms and the Fe plane is 1.41 
(1.40) A, and the Fe-Se-Fe angle for top (bottom) Se atoms is 108.5 
(108.7) degrees. The magnetic moment on each atom is calculated by 
integrating the spin density in Wigner-Seitz radius of 1.30 A, 1.16 A, 
2.14 A, 1.22 A, and 0.82 A for Fe, Se, Sr, Ti, and O atoms respectively. 
The calculated magnetic moment is 2.4 /n B on each Fe atom and is 
negligible on the other atoms. 

The energy bands for spin-up channel are shown in Fig. 2(a). The 
spin-down energy bands are the same as those for spin-up channel 
except for an exchange of X and X' points in BZ. There are two bands 
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Figure 2 | The spin-up energy bands of monolayer FeSe adsorbed on Ti0 2 
terminated SrTiO 3 (001) surface (a). The red and blue lines are contributed 
by 3d orbitals on Fe atoms and 2p orbitals on O atoms, respectively. The 
Fermi level is set to be zero. Panel (b) shows the influence of Hubbard- Uon 
the energy difference A r . 

[indicated as a and y in Fig. 2(a)] crossing the Fermi level E F near M 
point and foming an ideal pair of electron-hole bands with perfect 
Fermi-surface nesting, while around T point no metallic bands 
occur. Thus, the calculated Fermi surface only exists at the BZ corner, 
which agrees with the ARPES result in experiment 3 . In our calcula- 
tions, the antiferrodistortion of SrTi0 3 substrate does not have con- 
siderable influence to the electronic structure of adsorbed AFM FeSe 
monolayer. 

The disappearance of Fermi surface at BZ center is the result of a 
full electron filling of band fi. We find that the shape of band /? can be 
altered dramatically by changing the value of Hubbard- U correction. 
Actually, band /? shall cross the Fermi level at the BZ center in the 
absence of Hubbard- U correction. When the Hubbard- U correction 
is applied, band /? lowers its energy and a gap of 41 meV [indicated 
by green arrows in Fig. 2(a)] is obtained for U = 0.5 eV, which agrees 
with the experimental result 3 . The detailed relationship between the 
minimum energy difference of band ft and E F around the BZ center 
(denoted as A r ) and the Hubbard- U correction is shown in Fig. 2(b). 
One can see that A r and U are mainly in a linear relation. We have 
checked that this phenomenon does not exist in the paramagnetic 
FeSe/SrTi0 3 system, wherein there are energy bands crossing the 
Fermi level around the BZ center even for 5.0 eV Hubbard- U 
correction. 

Our PDOS (partial density of states) analysis shows that bands a 
and p are contributed by Fe 3d orbitals, while the hole-type band y 
comes from the O 2p orbitals. We analyze the wavefunctions of band 
P at r point and bands a, ft and y at M point. The wavefunction of 
band /? at T point shows a shape of a dumbbell in z- direction with a 
torus in xy-plane, which is the typical d z i feature. Whereas, the 
wavefunctions of bands a and ft at M point show vertical four-leaf 
shapes, which are the d xz and d yz features. The wavefunctions for 
band y are dumbbell shaped p x and p y orbitals that are mainly dis- 
tributed on the surface O atoms. The crossing between bands a and y 
indicates that charge transfer occurs from surface O atoms of SrTi0 3 
substrate to Fe atoms of FeSe monolayer. Our Bader analysis iden- 
tifies a charge transfer of -0.17 e per unit cell. 

Discussion 

From the above results, we know that the AFM FeSe monolayer is 
charge doped by the SrTi0 3 substrate. In order to analyze this effect, 
we study a free-standing AFM FeSe monolayer by doping charge 
from -0.5 e to 0.5 e. The results are shown in Fig. 3. We see that 
the shape of band /? changes dramatically when the system is posi- 
tively charged. Band p at T point crosses the Fermi level, and 
becomes a hole pocket. When the FeSe monolayer is negatively 
charged, the shape of band /? keeps static, while band a at T point 
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Figure 3 | The electronic energy bands of free-standing monolayer 
FeSe with 0.5 eV Hubbard- U correction and different levels of charge 
doping. The value of charge doping is changed from -0.5 eto 0.5 eper unit 
cell. The Fermi energy is set to zero. 

lowers its energy and finally crosses the Fermi energy, showing an 
electron pocket. Both positive and negative charge doping increase 
E F relative to the energy of band a and /? at M point. These results 
agree well with the previous study on FeSe monolayer without 
Hubbard- U correction 9 . However, as we have determined above, 
the realistic charge transfer in FeSe/SrTi0 3 is only -0.1 7e. By com- 
parison with Fig. 3, we find that this value is not large enough for the 
existence of electron pocket at T point. Thus, the energy bands of 
AFM FeSe monolayer are altered limitedly by pure charge doping 
from SrTi0 3 substrate. 

Besides the charge doping, another prominent feature in Fig. 2 is 
that bands a and ft are nonsymmetric along M-X and M-X' direc- 
tions, and a small splitting appears at M point. Whereas, for the free- 
standing AFM FeSe monolayer, bands a and /? are symmetric along 
M-X and M-X' directions, and they are degenerate at M point for 
both neutral and charged systems as shown in Fig. 3. Considering the 
net charge and dipole effect of SrTi0 3 substrate, a vertical electric 
field may exist in FeSe/SrTi0 3 . Then we further study the electric 
filed effect. In our further DFT calculations, a free-standing mono- 
layer FeSe is particularly designed to be exposed in a perpendicular 
electric filed, and the atomic positions are relaxed for each electric 
field strength. The calculation results are shown in Fig. 4. Clearly, 
when the electric field strength is zero, the energy bands are sym- 
metric along M-X and M-X' directions, and the bands of spin-up and 
spin-down components are the same. When the electric field 
strength is nonzero, however, we are astonished to see that the energy 
bands along M-X and M-X' directions become nonsymmetric, with 
an obvious spin splitting at M point. The splitting amplitude 
increases with increasing the electric field strength. So we verify that 
the energy bands near the Fermi energy and M point are distorted by 
the perpendicular electric field, and the energy of band a at X and X' 
points are different. Besides that, the features of spin-up and 
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Figure 4 | The electronic energy bands of free-standing monolayer FeSe 
with different vertical electric field strength. The value of electric field 
strength is changed from 0.0 eV/A to 1.0 eV/A. The energy bands near 
Fermi energy are distorted and show different behaviors for spin-up and 
spin-down channels. The Fermi energy is set to zero. 

spin-down bands are just opposite, namely, the energy bands in 
M-X' (M-X) direction for spin-up electrons are just the same as 
the energy bands in M-X (M-X') direction for spin-down electrons. 
Comparing the spin-up bands in Fig. 4 for free-standing FeSe mono- 
layer and those for FeSe monolayer on SrTi0 3 surface (Fig. 2), we see 
that the distortion of the energy bands are quite similar. Moreover, 
the FeSe monolayer energy bands near M point along T-M direction 
shows an energy gap of —30 meV when the applied electric field is 
0.5 eV/A, which agrees well with Fig. 3(i) and 3(j) in Ref.4. Thus, the 
effect of self- established electric field in FeSe/SrTi0 3 is confirmed. 

To explore the detailed mechanism of energy band distortion, we 
construct a tight-binding model Hamiltonian to describe the system. 
As we have discussed above, the corresponding wavefunctions of 
bands a and $ near M point are composed of d xz and d yz orbitals 
of Fe atom, and the wavefunction of energy band near T point 
mainly comes from d z i orbital of Fe atom. From Fig. 4 we know that 
the energy band distortion mainly occurs near M point. Thus, we will 
use d xz and d yz orbitals of Fe atoms in our tight-binding model, and 
neglect the d z i orbital for simplicity. However, the vertical electric 
field could not affect a model which only consists of d xz and d yz 
orbitals of Fe atoms, since all the centers of these orbitals are in a 
horizontal plane. Thus some out-plane orbitals should be included in 
our model. The most natural choice is to use the orbitals of top and 
bottom Se atoms. We assume the Se orbitals in our model are sym- 
metric in xy-plane, which may be Se 4p z orbital or its combination 
with Se 4s orbital. The following calculation shows that this assump- 
tion well describes the energy bands along M-X and M-X' directions. 

In our mean-field tight-binding model of AFM FeSe monolayer, 
we describe the spin-up and spin-down electrons separately. For 
spin-up electrons, the tight-binding model contains d xz and d yz orbi- 
tals of the spin-up Fe atom, besides one symmetric orbital for each Se 
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Figure 5 | (a) Atomic structure of FeSe monolayer. The red rectangle and green dash rectangle show two unit cells of AFM FeSe monolayer. Panel (b) 
shows the hopping schematic diagram of spin-up tight-binding model in a unit cell. The red, black, and gray circles show the positions of Fe atoms, top Se 
atoms and bottom Se atoms separately. In panel (b) the spin-down Fe atom is located at the cell center and is not plotted in this spin-up tight-binding 
model. The spin-down tight-binding model can be obtained by rotating panel (b) with 90 degrees or exchanging the positions of black and gray circles. 



atom. Similarly, for spin-down electrons, the tight-binding model 
contains d xz and d yz orbitals of the spin -down Fe atom, and one 
symmetric orbital for each Se atom. The hopping of spin-up elec- 
trons is schematically shown in Fig. 5. For spin-down electrons, 
whereas, the hopping schematic diagram can be obtained by switch- 
ing the positions of top and bottom Se atoms. 

Considering the symmetry of each orbital, the tight-binding model 
Hamiltonian for spin-up electrons is written as H^ up ^ = H^ ite + 

H^ p \ where H^ site and are shown in Eqs. (1) and (2) respect- 
ively. In these equations, c+ ? and c+ - are creation operators of Fe d xz 
and d yz orbitals, while c+ - and c£ ot - are those of top and bottom Se 
symmetric orbitals, respectively. For simplicity the spin indices are 
omitted. The parameter d is the distance 



top and bottom Se atoms are exchanged. Thus, the Hamiltonian for 
spin-down electrons is written as H^ down ^ =H^™te +H^ own \ where 

(down) nnA ^-(down) 



^on-site — ^ [ € d {^xzj C xz,r + C^j C yz ^j + (e 5 + dE)c^ p jC to pf 
r 

+ (£ s -dE)c£ ot jC bot ^ . 

r 

+ (cxzjCxzf+y + C yz,r C yz ,7 + x) 

+ fe (^xz,r C bot,r+x/2 + C yz,r C top,r+y/2^ ] + h.C. 

^on-site = Yl \^d( C xz,r C xzf + C yz,^ 
r 

+ (e s -dE)c+ t jC hot ^. 

^(down) _ I" / + , + \ 

— \j l [ C yz,r C yz,r+y^~ C xz,r C xz,f+xJ 
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H~ fe (^xz,r^ xz ^+y C yz,r C y z ^+^j 
+ fe (cxz,r C top,r+x/2 + ^z^botj+y/l^j ] +h.C. 



(2) 



(3) 



(4) 



between top (or bottom) Se atom and Fe plane, and E are the electric 
field strength. 

The tight-binding model Hamiltonian for spin-down electrons is 
similar to that of spin-up electrons, except for that the positions of 



^on-site' and H T 



are shown in equation (3) and (4) respectively. 
By performing the Fourier transformation, the model Hamilto- 
nians for spin-up and spin-down electrons in momentum space are 
given by Eqs. (5) and (6), 
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Figure 6 | Calculated tight-binding energy bands of monolayer FeSe 
for spin-up (a-c) and spin-down (d-f) channels. The energy bands 
distortion is shown near M point. The Fermi energy is set to zero. 
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where e dxz (k) =e d + t 1 cos(k y ) + t 2 cos(^), % z (k) = e d + ti cos(k x ) 
+ 1 2 cos (ky) , e top = e 5 + dE, and e bot = e 5 — dE. 

After diagonalization of Eqs. (5) and (6), the FeSe monolayer's 
energy bands near the Fermi surface are plotted in Fig. 6, in which 
we have used the parameters of e d = 0.36 eV, e 5 = — 4.57 eV, d = 1.37 
A, h = -0.01 eV, t 2 = 0.52 eV and t 3 = 0.71 eV. The results show the 
Fermi surface distortion clearly. The energy bands are symmetric 
along X-M-X', and are the same for spin-up and spin-down electrons. 
In the presence of a nonzero vertical electric field, an energy splitting 
appears at M point. The amplitude of this splitting increases with 
increasing the electric field strength. The energy bands become asym- 
metric along X-M-X'. These features are quantitatively consistent with 
our DFT results [see Fig. (4)]. Thus, our tight-binding model further 
explains the origin of energy bands distortion in FeSe/SrTi0 3 . Here, it 
should be noticed that in order to describe the electronic structures for 
the entire Brillouin zone, the tight-binding model should contains Fe 
d z 2 orbital besides the orbitals discussed above, since the energy band 
P in Fig. 2 at r point is mainly contributed by d z i orbital. 

To summarize, by performing the ab initio calculations, we have 
systematically investigated the electronic structure of AFM FeSe 
monolayer on Ti0 2 terminated SrTiO 3 (001) surface. We have shown 
that the interaction between SrTiO 3 (001) substrate and AFM FeSe 
monolayer induces charge transfer and a self- constructed electric 
field, which fundamentally modify the band structure and Fermi 
surface topology of the FeSe monolayer. In particular, the energy 
bands near the Fermi level are dramatically distorted by the electric 
field and show an obvious spin-resolved splitting around M point in 
BZ. As a result, the spin -resolved electron pockets evolved from band 
a become anisotropic and flat. We have proposed a tight-binding 
model Hamiltonian to reproduce these features. Also, we have shown 
that the present monolayer FeSe/SrTiO 3 (001) composite structure is 
an ideal electron-hole bilayer system, with electrons and holes 
respectively formed in FeSe monolayer and Ti0 2 . These findings 
are to be verified by the future spin-polarized ARPES experiments 
and are expected to shed light on the high-temperature supercon- 
ductivity or superfluidity at oxide interfaces. 

Methods 

The total energy and electronic structure calculations are performed by using projector 
augmented wave (PAW) method 15 " 18 . The exchange correlation potential is described 
by the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) 
type 19 . The plane wave cut-off energy is chosen to be 400 eV, which is converged in our 
test. 9X9X1 and 6X6X1 Monkhorst-Pack /c-points are used in the reciprocal space 
integration for 1 X 1 and \fl x \fl supercells. The system is relaxed until the force on 
each atom is smaller than 0.01 eV/A. To include the strong correlation effect for d 
electrons of Fe atoms, we adopt Hubbard- U correction. We use U = 0.5 eV for Fe-3d 
electrons in our main calculations, and the influence of different values of Hubbard- U 
on electronic structures are also discussed in the main text. All the DFT calculations are 
performed by using Vienna Ab-initio Simulation Package (VASP) 14 . 
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